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Ion  composition  provides  an  important  diagnostic  tool  for  identifying  and 
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magnetospheric  charged  particles.  However,  the  situation  is  not  as  simple  as 
had  been  envisioned  a  decade  ago.  Ring-current  and  radiation-belt  ions  are 
believed  to  come  from  either  the  sun  or  the  earth's  ionosphere  or  both.  Entry 
is  not  direct  in  either  case.  Solar-wind  velocity  and  ion  composition  vary 
considerably  with  time.  Ion-velocity  distributions  are  partially  thermalized 
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-in  the  magnetosheath  after  crossing  the  bow  shock.  Entry  of  solar-wind  ions 
into  the  magnetosphere  is  modulated  in  part  by  reconnection  efficiency,  i.e., 
by  the  ratio  of  cross-tail  electric  field  to  "asymptotic”  interplanetary  elec¬ 
tric  field.  Non-adiabatic  motion  in  the  plasma  sheet  may  enable  the  cross¬ 
tail  electric  field  to  heat  the  various  ion  species  by  different  amounts. 

Entry  of  ionospheric  ions  into  the  plasma  sheet  seems  to  occur  primarily 
through  afternoon-  and  evening-sector  auroral  arcs,  rather  than  through  the 
polar  wind  as  had  previously  been  postulated.  Energization  occurs  through 
Interaction  with  the  auroral  potential  structure  (which  varies  with  magnetic 
activity)  and  through  plasma-sheet  processes  cited  above. ^=4hes£__are  some  of 
the  uncertainties  that  affect  the  outer  boundary  condition  for  the  pfobl^m^of 
magnetospheric  ion  composition.  Radial  transport  of  magnetospheric  ions  at 
ring-current  energies  (E  ~  10-100  keV)  is  achieved  primarily  by  unsteady  elec¬ 
trostatic  convection.  The  corresponding  diffusion  coefficient  is  independent 
of  charge  and  mass  among  particles  having  drift  periods  »  2n  times  the  postu¬ 
lated  decay  time  of  an  impulse  in  the  convection  electric  field.  Erosion  of 
the  ring  current  occurs  primarily  through  charge  exchange  and  perhaps  through 
wave-particle  interactions.  Charge  exchange  favors  the  survival  of  0+  and  He+ 
over  H+  and  He-14"  at  energies  <  50  keV,  but  account  must  be  taken  of  population 
disparities  in  the  source  region  and  of  time  scales  associated  with  radial 
diffusion  therefrom.  Colllsional  processes  are  less  Important  at  radiation 
belt  energies  (E  >  200  keV),  and  the  mapping  of  ion  distributions  from  the 
source  region  to  low  L  values  can  be  partially  understood  in  terms  of  a  common 
diffusion  profile  with  M  and  J  (first  two  adiabatic  invariants)  conserved. 

The  pitch-angle  anisotropy  of  ring-current  ions  is  likely  to  generate  electro¬ 
static  ion-cyclotron  waves  outside  the  plasmasphere  and  electromagnetic  ion 
cyclotron  waves  inside  the  plasmasphere.  These  should  result  in  preferential 
pitch-angle  diffusion  of  the  major  ionic  constituent  and  in  a  possibly  weaker 
(parasitic)  pitch-angle  diffusion  of  the  minor  ionic  constituents. 
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FIGURES 


Radial  variation  of  driven  (a)  by  step-like  magnetic 
impulses  and  (b)  by  electrostatic  impulses  that  rise  quickly 
and  decay  slowly  on  the  time  scale  of  an  ion's  drift  period... 

Empirical  ^-folding  energies  of  observed  proton  spectra 
(Davis  and  Williamson,  1963)  at  equatorial  pitch  angles 
consistent  with  (5)  and  mapping  to  aQ  ■  10°,  20°,  30°,  and 
90°  (respectively)  at  L  ■  7;  expected  variation  of  energy 
with  L  for  individual  protons  having  constant  M  and  J,  for 
selected  values  of  energy  and  equatorial  pitch  angle  at 
L  -  7  (Dungey  et  al.,  1963) . 

Equatorial  pitch-angle  distributions  of  alpha  particles  and 
protons  in  the  same  four  energy/nucleon  passbands,  as  observed 
on  spacecraft  0V1-19  (Blake  et  al.,  1973) . 

Charge-exchange  lifetimes  against  neutralization  in  model 
atomic-hydrogen  atmosphere  for  equatorlally  mirroring  (K  -  0) 
protons  and  helium  ions  at  selected  values  of  M/A,  measured 
in  MeV/gauss-nucleon  (Cornwall,  1971;  Schulz  and  Lanzerotti, 
1974) 

Products  of  ion  velocity  and  cross  section  for  charge  exchange 
in  atmosphere  of  atomic  hydrogen  (Cornwall  and  Schulz,  1979), 
based  on  compilations  summarized  by  Cornwall  (1971)  and 
Tinsley  (1976) . 


INTRODUCTION 


In  more  innocent  times  it  was  believed  that  the  ion  composition  of  the 
ring  current  and  radiation  belts  at  a  specified  energy/nucleon  should  be  the 
same  as  the  ion  composition  of  the  solar  wind.  The  idea  was  that  traversal  of 
the  bow  shock  would  separately  thermalize  each  ionic  constituent  of  the  solar 
wind  so  as  to  produce  a  Maxwellian  distribution  having  a  temperature  ~  1 
keV/nucleon  as  an  outer  boundary  condition  at  L  ~  10  on  the  phase-space  densi¬ 
ty,  and  that  the  diffusion  equation  would  impose  a  common  profile  (defined  by 
the  requirement  of  a  divergence- free  radial-diffusion  current)  on  all  the 
constituents  in  the  ring  current  and  radiation  belts.  This  concept  has  large¬ 
ly  failed  to  be  supported  by  the  observational  data  and  so  is  considered 
obsolete. 

One  difficulty  with  the  above  concept  is  that  the  enfolding  energy  pro¬ 
vided  by  thermalizlng  the  solar  wind  is  too  smal*  by  an  order  of  magnitude 
(see  below)  for  assignment  to  a  boundary  at  L  ~  10.  Another  difficulty  is  the 
neglect  of  ionospheric  ions,  which  can  enter  the  magnetosphere  with  several 
keV  of  energy  by  virtue  of  the  field-aligned  potential  drop  characteristic  of 
auroral  arcs  in  the  afternoon  and  evening  sectors.  It  seems,  however,  that 
entry  into  the  radiation  belt  or  ring  current  is  not  direct  in  either  case, 
but  rather  occurs  via  the  plasma  sheet.  Violation  of  the  first  adiabatic 
invariant  of  charged-particle  motion  within  the  plasma  sheet  may  enable  the 
cross-tail  electric  field  to  heat  the  various  ionic  species  further,  but  not 
necessarily  by  the  3ame  amount  or  by  the  same  amount  per  nucleon.  Moreover, 
entry  of  solar  wind  ions  into  the  plasma  sheet  is  modulated  in  part  by  recon¬ 
nection  efficiency,  i.e.,  by  the  ratio  of  cross-tail  electric  field  to  asymp¬ 
totic  interplanetary  electric  field.  The  reconnection  efficiency  can  vary 
considerably  with  time,  as  can  the  solar-wind  velocity  and  ion-composition. 


Likewise,  the  field-aligned  auroral  potential  drop  can  vary  considerably  with 
time.  These  are  some  of  the  uncertainties  that  can  affect  the  outer  boundary 
condition  for  the  problem  of  magnetospheric  ion  composition  as  a  function  of 
energy/nucleon. 

Radial  transport  of  magnetospheric  ions  at  ring-current  energies  (C  ~  10 

-100  keV)  is  achieved  primarily  by  unsteady  electrostatic  convection.  The 

corresponding  diffusion  coefficient  is  Independent  of  charge  and  mass  among 

twenty  minutes  is 

particles  having  drift  periods  »  2  hr  if  the  postulated  decay  time  of  an 
impulse  in  the  convection  electric  field.  However,  the  diffusion  current  is 
not  divergence-free,  since  distributed  losses  constitute  an  essential  ingredi¬ 
ent  of  ring-current  dynamics.  Erosion  of  the  ring  current  occurs  primarily 
through  charge  exchange  and  perhaps  through  wave-particle  interactions. 
Charge  exchange  favors  the  survival  of  0+  and  He+  over  H+  and  He++  at  ener¬ 
gies  5  50  keV,  but  account  must  be  taken  of  population  disparities  in  the 
plasma  sheet  and  of  time  scales  associated  with  radial  diffusion  therefrom. 
Thus,  it  is  possible  for  H+  to  predominate  over  0+  and  He+  in  the  early  stages 
of  ring-current  formation  but  not  in  the  later  stages  of  a  magnetic  storm. 

Collislonal  processes  are  less  Important  at  radiation-belt  energies  (E 
200  keV),  and  the  mapping  of  ion  distributions  from  the  outer  boundary  to  low 
L  values  can  be  partially  understood  in  terms  of  a  common  diffusion  profile 
with  tl  and  J  (first  two  adiabatic  invariants)  conserved.  One  consequence  of 
this  is  the  appearance  of  very  large  ionic  pitch-angle  anisotropies  in  the 
inner  radiation  zone.  Indeed,  the  pitch-angle  anisotropy  at  low  L  values  is  a 
reflection  of  the  energy  spectrum  at  the  outer-boundary  location.  Inner-zone 
observations  typically  reveal  that  the  alpha-particle  anisotropy  is  even 
larger  than  the  proton  anisotropy  at  the  same  energy /nucleon.  This  means  that 
the  a/p  (alpha/proton)  Intensity  ratio  must  vary  inversely  with  magnetic 


latitude  along  a  field  line,  as  indeed  it  does.  The  occurrence  of  a  larger 
alpha-particle  chan  proton  anisotropy  at  fixed  energy/nucleon  can  be 
understood  in  terms  of  an  alpha-particle  temperature  chat  is  less  chan  four 
times  the  proton  temperature  at  Che  outer  boundary  of  crapped  radiation. 

Much  of  the  observational  research  emphasis  on  ion  composition  in  recent 
years  has  focused  on  energies  S  20  keV,  i.e.,  on  energies  so  low  that  they  are 
not  truly  in  the  radiation- belt  range.  Indeed,  such  energies  are  so  low  that 
they  hardly  include  a  significant  portion  of  even  the  ring-current  spectrum. 
It  is  therefore  quite  important  to  avoid  sweeping  generalizations  in  the 
interpretation  of  ion-composition  results.  Meanwhile,  progress  (e.g., 
Gloeckler,  1977)  is  being  made  toward  the  construction  of  ion-identifying 
spectrometers  that  operate  in  what  had  been  the  energy  "gap"  between  ~  20 
keV/ charge  and  ~  100  keV/nucleon,  as  required  for  proper  scrutiny  of  the  ring- 
current  spectrum. 

The  purpose  of  this  work  is  to  outline  the  basic  principles  upon  which 
•  the  theoretical  study  of  magnetospheric  ion  composition  is  founded.  The 
author  helped  to  compile  a  much  more  thorough  review  of  the  literature  on  this 
subject  only  a  few  years  ago  (Cornwall  and  Schulz,  1979).  Not  much  would  be 
gained  by  providing  the  results  of  yet  another  literature  search  at  this 
time.  Thus,  the  focus  of  the  present  review  is  directed  instead  on  the  endur¬ 
ing  fundamentals.  These  relate  to  the  problems  of  access,  energization, 
transport,  and  loss  of  protons  and  other  ions  with  respect  to  the  earth's 
magnetosphere  and  radiation  belts. 


ACCESS 


The  main  sources  of  magnecospheric  ions  are  the  galaxy,  the  sun,  and  the 
ionosphere.  The  galaxy  provides  cosmic-ray  ions  that  are  themselves  too 
energetic  to  be  trapped  by  the  earth's  magnetic  field.  However,  the  backscat- 
ter  (albedo)  products  of  spallation  reactions  that  occur  when  such  cosmic-ray 
ions  strike  the  upper  atmosphere  include  many  neutrons  in  the  10-100  MeV 
range.  The  occasional  decay  of  a  neutron  from  this  population  contributes  a 
proton  in  the  10-100  MeV  energy  range  and  an  electron  of  energy  <  800  keV  to 
the  radiation  belts.  Other  sources  of  electrons  are  much  more  impor  than 
this  one.  Other  sources  of  inner-zone  ions  in  the  10-100  MeV  energy  range  are 
negligible.  Thus,  inner-belt  ion  composition  is  not  much  of  a  problem  at 
these  very  high  energies:  the  ions  in  this  range  are  almost  exlusively  pro¬ 
tons.  The  source  described  in  this  paragraph  is  known  by  the  acronym  CRA1JD 
(cosmic-ray-albedo-neutron  decay). 

Drift  shells  at  synchronous  altitude  and  beyond  are  occasionally  popu¬ 
lated  by  solar-flare  ions  (presumably  protons)  of  energy  5-70  MeV  (Paulikas 
and  Blake,  19b9).  The  mode  of  access  in  this  case  is  largely  direct,  in  that 
the  drift  shells  populated  are  chose  that  pass  within  a  gyro-radius  of  the 
magnetopause  or  neutral  sheet. 

Of  greater  interest  in  the  context  of  magnetospheric  ion  composition  are 
particles  of  much  lower  energy,  perhaps  E  ~  1-10  keV,  for  which  the  modes  of 
access  are  largely  indirect.  These  include  magnetosheach  ions  derived  from 
the  shocked  solar  wind  and  auroral  ions  energized  by  the  electric  potential 
drop  that  occurs  along  magnetic  field  lines  in  the  PM  (afternoon-evening) 
sector  of  the  auroral  oval.  Magnetosheath  ions  with  E  ~  1  keV  might  impinge 
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on  the  magnetosphere  anywhere,  but  the  convection  electric  field  tends  to 

expel  them  from  the  dayside  region  of  closed  field  lines.  Conversely,  charged 

particles  thac  enter  at  the  polar  cleft  (a  dayside  region  of  wea*.  magnetic 

field  surrounding  the  neutral  points)  are  drawn  into  the  tail  lobes  by  the 
convection  electric  field,  as  are  particles  that 

impinge  on  the  tail-lobe  surfaces.  Characteristics  of  the  resulting  plasma 

mantle  have  been  described  by  Sckopke  and  Paschmann  (1978). 

The  convection  velocity  of  plasma  in  either  tail  lobe  is  smaller  than  the 
solar-wind  velocity  by  a  factor  of  order  eB  /B  «  1,  where  Bt  is  the  tail 
field,  Bz  is  the  southward  component  of  the  interplanetary  magnetic  field,  and 
e  is  the  "reconnection  efficiency"  defined  by  Kennel  and  Coroniti  (1975)  as 
the  ratio  of  the  east-west  component  of  the  cross-trail  electric  field  to  that 
of  the  interplanetary  electric  field.  Observations  of  Burton  et  al.  (1975) 
suggest  a  "half-wave  rectifier*  model  such  that  e  is  practically  zero  for 
northward  Bz  but  at  least  0.2  (and  perhaps  even  -  1)  for  southward  Bz>  The 
factor  cB  / B  is  quite  small  (S  0.1)  in  any  case,  and  so  only  the  extreme 

Z  L 

flanks  of  the  plasma  mantle  (i.e.,  regions  already  within  a  few  earth  radii  of 
the  current  sheet  in  the  tail)  are  likely  to  populate  the  plasma  sheet  by 
convection  within  the  first  few  tens  of  earth  radii  downstream  from  the 
earth.  The  reversal  of  magnetic-field  direction  across  the  current  sheet  and 
the  dawn-to-dusk  direction  o:  the  cross-tail  electric  field  make  the  low- 
latitude  A11  (pre-dawn)  flank  an  inviting  place  for  magnetosheath  ions  to  enter 
the  oagnetosphere  in  any  case. 

The  access  of  ionospheric  ions  to  the  plasma  sheet  is  made  possible  by 
the  electrostatic  potential  structure  characteristic  of  the  PM  (afternoon- 
evening)  sector  of  the  auroral  oval.  It  is  in  this  region  (near  the  boundary 
between  closed  and  open  magnetic  field  lines)  that  there  develops  a  virtual 
discontinuity  in  the  magnetospheric  convection  electric  field  if  one  treats 


magnetic  field  lines  as  electrostatic  equipotentials.  There  is  a  Pedersen 
conductance  owing  to  the  collisionality  of  the  ionosphere  and  an  Ohmic  impe¬ 
dance  due  to  the  mirror  forces  exerted  by  the  inhomogeneous  magnetic  field  on 
the  hoc  auroral  plasma.  It  seems  that  these  properties  conspire  with  the 
requirement  of  current  conservation  (Ampere's  law)  to  produce  a  multi-kilovolt 
potential  drop  along  the  field  lines  of  the  auroral  oval.  The  geometry  is 
such  that  the  corresponding  parallel  (to  B)  component  of  the  electric  field 
should  be  upward  in  the  PM  sector  and  downward  (if  present)  in  the  All 
sector.  Statistics  compiled  by  Ghielmetti  et  al.  (1978)  from  S3-3  data  sug¬ 
gest  that  upgoing  ion  beams  are  about  six  times  more  probable  in  the  PM  sector 
chan  in  the  AM  sector,  and  about  44  times  more  probable  in  the  dusk  quadrant 
(15-21  hr,  magnetic  local  time)  than  in  the  dawn  quadrant  (03-09  hr,  MLT) . 

The  contribution  of  auroral  ion  beams  to  the  plasma  sheet  presumably 
exceeds  that  of  the  polar  wind,  which  had  been  until  recently  the  favored 
medium  (e.g.,  Axford,  1970)  for  supplying  ionospheric  ions  to  the  plasma 
sheet.  The  polar  wind  theoretically  consists  of  protons  (along  with  a  few 
helium  ions)  having  energies  £  10  eV.  The  auroral  ion  beams  seem  to  be  rich 
in  oxygen  as  well,  with  ion  energies  of  several  keV  at  least.  It  is  true  chat 
the  polar  wind  delivers  ions  from  the  entire  polar  cap  rather  than  from  only 
the  PM  half  of  the  auroral  oval,  but  the  deficiencies  in  ion  energies  and 
oxygen  ions  now  seem  to  make  the  polar  wind  a  less  important  medium  (tnan  the 
auroral  ion  beams)  for  the  transfer  of  ionospheric  ions  to  the  plasma  sheet. 
Thus,  it  seems  that  the  main  sources  of  1-10  keV  ions  for  the  plasma  sneet  are 
(1)  magnetosheath  plasma  derived  from  the  shocked  solar  wind  and  (2)  iono¬ 
spheric  plasma  from  the  PM  sector  of  the  auroral  oval. 
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energization 


The  importance  of  Che  plasma  sheet  for  magnetospheric  ion  composition 
resides  in  the  adiabatic  and  non-adiabatic  processes  by  which  ions  are  further 
energized  there  before  being  transported  into  the  region  of  closed  drift 
shells.  The  adiabatic  energization  process  entails  the  usual  gradient- 
curvature  drift  of  nightside  ions  across  equipotentials  of  the  dawn-to-dusk. 
electric  field,  or  (equivalently)  the  sunward  convection  of  nightside  ions 
into  a  spatially  increasing  magnetic-mirror  field.  This  process  occurs  in  the 
earthward  portion  of  the  nightside  plasma  sheet,  i.e.,  in  a  region  where  the 
magnetic  field  lines  are  topologically  dipolar  but  from  which  the  adiabatic 
drift  shells  intersect  the  magnetopause.  The  particles  in  question  are  "quasi- 
trapped"  in  the  terminology  of  Roederer  (1970),  but  the  failure  of  their  drift 
shells  to  close  is  as  much  the  fault  of  the  convection  electric  field  as  it  is 
of  the  gradient-curvature  drift. 

The  non-adiabatic  energization  process  occurs  in  the  tailward  portion  of 
the  plasma  sheet,  where  the  northward  (z)  component  of  the  earth's  magnetic 
field  is  too  weak  to  enforce  conservation  of  the  first  two  adiabatic  invari¬ 
ants.  A  plasma  physicist's  prototype  for  the  magnetic  field  in  this  distinct¬ 
ly  non-dipolar  region  might  be  expressed  as 

A  A 

B  ■  z  B  -  x  B  tanh  (z/z_),  (1) 

"A  Z  TL—lTl  L  U 

where  x  is  the  anti-sunward  unit  vector,  Zq  is  the  “half-width"  of  tne  current 
sheet,  and  »  Bz«  A  model  electric  field  consistent  with  (1)  is  the  uniform 
field  E  ■  -  y|Ey|.  A  theorem  of  Stern  and  Palmadesso  (1975)  asserts  that 
there  is  no  net  gradient-curvature  drift  in  a  magnetic  field  such  as  (1)  if 
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B2,  Bt,  and  Zq  are  constants.  The  proof  follows  from  the  construction  of  a 
bounce-averaged  Hamiltonian  function  that  depends  on  the  first  two  adiabatic 
invariants  M  and  J  but  not  on  the  coordinates  xq  and  yQ  at  which  the  guiding 

field  line  crosses  the  plane  z  ■  0.  An  immediate  corollary  is  that  the  elec- 

a  ^ 

trie  field  ■  -  yiE^j  produces  a  bounce-averaged  drift  such  that  Xq  » 
-(c/Bz))Eyl,  but  that  particle  energization  does  not  occur. 

The  theorem  of  Stern  and  Palmadesso  (1975)  clearly  does  not  apply  if  Bz  - 

* 

0.  In  this  case  there  is  a  gradient  drift  in  the  +  y  direction  for  ions 

adiabatic  with  respect  to  M,  but  no  curvature  drift  since  the  field  lines  are 

* 

rectilinear.  There  is  also  a  net  convection  in  the  z  direction  (toward  the 

neutral  sheet)  for  all  ions  and  a  net  drift  in  the  -  y  direction  for  ions  that 

MMA 

fail  to  conserve  M  because  their  trajectories  cross  the  plane  z  *  0.  More¬ 
over,  the  ions  that  thus  fail  to  conserve  M  are  energized  by  the  uniform 
electric  field.  The  theorem  fails  in  the  case  Bz  •  0  because  the  second 

adiabatic  invariant  J  no  longer  exists,  or  (equivalently)  because  the  bounce 

average  is  no  longer  defined,  even  for  particles  that  conserve  the  first 
invariant  M. 

The  theorem  of  Stern  and  Palmadesso  (1975)  likewise  fails  even  for  B„  *  0 

«* 

if  the  particle  of  interest  has  too  large  a  rigidity  to  conserve  M  and  J  in 
the  field*  configuration  given  by  (1).  Thus,  the  dawn-to-dusk  electric  field 
can  be  expected  to  add  extra  energy  to  ions  that  have  entered  the  plasma  sheet 
with  rigidity  p/q  >  ZqBz/c.  Recent  results  of  Lyons  and  Speiser  (1982) 

suggest  that  the  energy  gained  in  this  way  is  Indeed  an  increasing  function  of 
initial  ion  energy.  This  is  a  geophysical  analogy  of  the  rich  becoming 
richer.  Moreover,  since  the  energization  mechanism  involves  only  single- 
particle  motion  in  crossed  electric  and  magnetic  fields,  tnere  is  no  expecta¬ 
tion  that  any  of  the  energy  gained  by  the  favored  ions  of  high  initial  energy 


will  "trickle  down"  to  the  less  favored  ions  of  low  initial  energy. 

The  above-described  idea  for  charged-particle  energization  in  the  tail- 
ward  portion  of  the  plasma  sheet  was  first  worked  out  in  quantitative  detail 
by  Speiser  (1965).  It  is  quite  remarkable  that  such  a  perceptive  work  has 
received  so  little  attention  over  the  past  seventeen  years.  In  any  event,  the 
occurrence  of  an  "ion-heating"  process  based  on  single-particle  motion  makes 
it  premature  to  invoke  plasma  turbulence  or  reconnection  as  the  energization 
mechanism  of  the  plasma  sheet.  It  appears  from  the  work  of  Lyons  and  Speiser 
(1982)  that  non-adiabatic  single-particle  motion  in  the  crossed  electric  and 
magnetic  fields  is  sufficient  to  account  for  ion  temperatures  observed  in  the 
tailward  portion  of  the  plasma  sheet.  Still  open  is  the  important  question  of 
whether  the  various  ionic  species  are  "heated"  by  different  amounts  (per 
particle,  per  unit  charge,  or  per  nucleon)  by  this  mechanism.  Sunward  convec¬ 
tion  by  the  dawn-to-dusk  electric  field  delivers  the  ions  thus  energized  into 
the  earthward  portion  of  the  plasma  sheet,  where  further  energization  occurs 
through  the  conservation  of  M  and  J  in  the  crossed  Inhomogeneous  magnetic  and 
electric  fields  that  prevail  there.  Of  course,  a  charged  particle  can  gain 
even  more  energy  if  it  is  transported  (e.g.,  by  guiding-center  diffusion) 
across  the  boundary  between  open  (plasma-sheet)  and  closed  (ring-current) 
drift  shells,  and  from  there  to  progressively  smaller  L  values  within  the 
region  of  geomagnet ically  trapped  particles.  However,  the  subject  of  charged- 
particle  energization  by  guiding-center  diffusion  is  Included  below  in  the 
section  about  charged-particle  transport  across  adiabatic  drift  shells. 


TRANSPORT 


For  charged  particles  that  conserve  their  first  two  adiabatic  in¬ 
variants  M  and  J,  the  problem  of  non- adiabatic  transport  typically 
entails  violation  of  the  third  invariant  #  through  radial  diffusion.  At 
radiation-belt  energies  it  is  usual  to  adopt  the  dimensionless  label  L 
(roughly  speaking,  the  equatorial  radius  of  the  drift  shell  in  earth 
radii)  as  diffusion  coordinate,  and  this  can  be  done  by  means  of  the 
transformation  JL  *  2ir(p/a)i  ®|"1-,  where  p  is  the  magnetic  moment  of 
the  earth  and  a^  is  the  radius  of  the  earth.  The  diffusion  equation  for 
the  drift-averaged  phase-space  density  T  then  takes  on  the  form 


(2) 


where  diffusion  coefficient  (e.  g.  ,  Schulz  and  Lanzerotti, 

1974).  The  functional  form  of  is  necessarily  model-dependent, 

but  the  usual  treatment  yields  results  of  the  form  +  ^LL 


where 


D(.n?)  «  7  X  10'9L10[Q(y)/180D(y)]2  day'1 


represents  the  contribution  from  step-like  impulses  in  the  magneto - 
spheric  B  field  and 

as  1  X  10*10L10(yZM0y2/M)2[T(y)/2D(y)]2[  1  +  (Ajt)'2]'1  day'1 

(3b  1 

represents  the  contribution  from  exponentially  decaying  impulses  in 
the  electrostatic  field  that  drives  magnetospheric  convection.  The 
auxiliary  functions  Q(y).  D(y),  and  T(y)  are  well  approximated  (Schulz 


and  Lanzerotti,  1974;  Davidson,  1976)  by  the  algebraic  forms 


Q(y)  as  -  27.12667  -  45.  39913y4  +  5.88256y8  (4a) 

D(y)  «  0.4600577  +  0.  1066 1 54  y3 /4  -  0.  1997662y  (4b) 

and 

T(y)  a:!  1.  3801730  -  0. 6396925  y3  /4  (4c) 


where  y  (=  sin  a^)  is  the  sine  of  the  equatorial  pitch  angle.  Other 
symbols  appearing  in  (3)  include  the  ratio  y  relativistic  mass  to 
rest  mass,  the  ionic  charge-state  number  Z,  the  scale  factor  Mq  »  1 
GeV/gauss,  the  ionic  drift  frequency  and  the  decay  time  t  (  — 

1200  sec)  of  a  model  electrostatic  impulse.  The  factor  [1  +  (ft^i")"2] 
differs  by  at  most  2.5%  from  unity  for  ions  having  drift  periods  2ir/£l, 
S  20  min.  i.  e. ,  for  ions  of  kinetic  energy  E  2  1.2Z/L  MeV,  and  so 
the  factor  [1  +  (D^t)-2]  in  (3b)  may  be  regarded  as  essentially  con¬ 
stant  for  such  particles. 


Figures  la  and  lb  show  the  variation  of  and  with  L 


in  the  limit  S2^T  >>  1  for  nonr elativistic  ions  having  the  same  energy 
(and  selected  values  of  y)  at  L  =  7.  The  reader  may  be  surprised  to 
see  that  scales  as  L3*3  only  for  y  *  1.  For  Sq  *  90°  account 

must  be  taken  of  the  variation  of  y  with  L  at  fixed  M  and  J.  Accord¬ 
ing  to  Chen  and  Stern  (1975)  this  variation  is  well  approximated  by 
the  algebraic  form 


where 


1  +  1.  35048  X  -  0.030425  X4/3 

+  0.  10066  X5 /3  +  (X/2.  760346)2 

(La/fi)(J3/8mQM)  for  a  particle  of  rest  mass  m». 


(5) 


The 


Figure  1.  Radial  variation  of  driven  (a)  by  step-like  magnetic 

impulses  and  (b)  by  electrostatic  impulses  that  rise  quickly  and  decay 
slowly  on  the  time  scale  of  an  ion's  drift  period.  Impulses  are  pre¬ 
sumed  to  conserve  M  and  J  so  that  y  varies  in  accordance  with  (5)  as 
a  function  of  L  for  each  value  of  y^  (denoting  the  value  of  y  at  L  *  7), 
i.  e,  for  y^  *  0,  0.4,  0.6,  and  1.0  in  (a)  and  for  y^  =  0,  0.2,  0.4, 
0.6,  0.8,  and  1.0  in  (b).  Dashed  lines  (y^  =  0)  are  not  realized  in 
practice,  because  of  the  loss  cone  (Schulz  and  Lanzerotti,  1974).  ' 


limiting  case  y  =  0  (X  *  oc)  thus  yields  y~  oc  1  /L  in  (3b),  from  which 
it  follows  that  scales  only  as  L3  in  this  limit  (see  Figure  lb, 

dashed  line).  On  the  other  hand,  it  seems  from  Figure  la  that  the 
deviation  of  from  an  scaling  is  almost  negligible.  Thus, 

for  ion  energies  sufficiently  high  that  charge  exchange  and  Coulomb 
drag  can  truly  be  neglected  as  they  are  in  writing  (2),  one  might  an¬ 
ticipate  a  steady-state  diffusion  profile  of  the  form 


7(M,  J,  L) 


'  1  -  (l^/L)7  ' 

*  y 

1  -  (Lj/L  )' 


7(M,  J,  L*) 


(6) 


j*e 

where  L  (—  10)  represents  the  outermost  closed  drift  shell  and  Lj 
represents  the  innermost  drift  shell  that  can  trap  particles  of  a  given 
K  *  J  /Sm^M.  Thus,  the  entire  energy  spectrum  of  7(M,  J,  L)  at  L 
<  L  is  determined  by  the  boundary-value  spectrum  contained  in  7(M, 
J,  L  ).  Of  particular  interest  in  the  case  of  nonrelativistic  ions  is 
the  specification 

7(M,  J,  L*)  =  [  (y*)2(L*a)3E*/)iM]i+1  7*  exp[  -  (AM/(y*)2(L*a)3E* ] 

(7) 

where  1,7,  and  Eq  are  constants.  The  special  case  l  ~  0  corre¬ 
sponds  to  an  exponential  spectrum,  since  the  particle  flux  (per  unit 

energy  and  solid  angle)  is  equal  to  2m^E?.  The  substitution  piM  = 

3  3  2 

L  a  y  E  in  (7)  is  valid  for  all  L  values  in  the  dipolar  field  model 
and  yields 


7(M,  J,  L*)  =  [(y*/y)2(L*/L)3(E* /E)]/  +  1  7* 

X  exp[-  (y/y*)2(L/L*)3(E/E*)]  (8) 


which  can  be  written  more  compactly  as 


(9a) 


KM.J.L*)  =  (Eq/E)1^1  T  exp(  -  E/Eq) 

by  introducing  the  parameter 

Eq  ■  (y*/y)2(L*/L)3E*  (9b) 

It  follows  from  (9)  that  the  index  l  of  a  power-law  spectrum  remains 
invariant  through  the  diffusion  profile  given  by  (6),  whereas  the  e- 

>fs 

folding  energy  Eq-  of  an  exponential  spectrum  varies  with  L  and  y  in 
the  same  manner  as  the  energy  of  an  individual  particle  for  which  M 
and  J  are  conserved-  Moreover,  the  variation  of  Eq  with  L  is 
stronger  for  y  —  1  than  for  y  «  1,  and  this  property  creates  an¬ 
isotropy  in  the  pitch-angle  distribution  observed  at  fixed  E,  over  and 

7  ♦  7 

above  the  anisotropy  contained  in  the  factor  [1  -  (L^/L)  ] /[  1  -  (L^/L  )  ], 
at  L  <  L*. 

Application  of  the  foregoing  ideas  to  a  body  of  observational  data 
is  illustrated  in  Figure  2  (Dungey  et  al.  ,  1965).  The  solid  curves 
were  constructed  by  measuring  the  £- folding  energies  Eq  of  exponen¬ 
tial  proton  spectra  (Davis  and  Williamson,  1963)  at  equatorial  pitch 
angles  corresponding  [via  (5)]  to  o-q  =  10°,  20°,  30°,  and  90°  at  L  * 

7.  The  dashed  curves  show  the  corresponding  energy  variation  of  an 
individual  proton  but  are  positioned  on  the  logarithmic  ordinate  so  as 
to  best  fit  the  observed  variation  of  Eq  with  L.  The  dashed  curves 
seem  to  converge  on  a  common  (isotropic)  value  of  Eq  **•  12  keV  at 

ajs 

L  **  10.  It  is  natural,  therefore,  to  associate  the  parameters  Eq  ~ 

12  keV  and  L  —  10  with  the  outer  boundary  condition  on  the  proton 
radiation  belt,  on  which  the  observations  of  Davis  and  Williamson 
(1963)  spanned  the  energy  range  E  *  0.1-1.  6  MeV.  These  results  for 
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Figure  2.  Empirical  enfolding  energies  (solid  curves)  of  observed  pro 
ton  spectra  (Davis  and  Williamson,  1963)  at  equatorial  pitch  angles 
consistent  with  (5)  and  mapping  to  aQ  =  10°,  20°,  30°,  and  90°  (re¬ 
spectively)  at  L  =  7;  expected  variation  (dashed  curves)  of  energy  with 
L  for  individual  protons  having  constant  M  and  J,  for  selected  values 
of  energy  and  .equatorial  pitch  angle  at  L  =  7  (Dungey  et  al. ,  1965). 


£q  and  L  suggest  that  the  adiabatic  and  non- adiabatic  energization 
processes  operative  in  the  plasma  sheet  must  be  able  to  provide  a 
proton  temperature  —  12  keV  at  L  —  10,  at  least  for  the  high-energy 
tail  of  the  proton  distribution  there.  The  appearance  of  a  boundary 
at  L  —  10  between  closed  and  open  drift  shells  is  quite  plausible  for 
the  (high)  proton  energies  of  interest. 

The  large  anisotropy  inherent  in  (8)  and  (9)  for  L  «  L*  offers 
a  likely  explanation  for  the  occurrence  of  pitch-angle  distributions 
such  as  those  illustrated  in  Figure  3  (Blake  et  al.  ,  1973).  The  pro¬ 
ton  anisotropy  at  L  =  1.  85  is  an  increasing  function  of  energy,  as 
would  be  expected  from  the  presence  of  an  exponential  factor  in  (9). 

The  alpha-particle  anisotropy  is  even  larger  than  the  proton  anisotropy 
at  the  same  energy /nucleon,  as  would  be  expected  if  the  alpha-particle 
temperature  were  somewhat  less  than  four  times  the  proton  tempera- 
ture  at  the  outer  boundary  of  trapped  radiation,  i.  e.  ,  at  L  =  L  . 

Such  a  boundary  condition  is  not  really  implausible  in  view  of  the 
adiabatic  and  non-adiabatic  plasma- sheet  processes  (see  above)  that 
must  intervene  before  a  solar -wind  ion  enters  the  radiation  belts. 
Moreover,  there  is  no  expectation  that  plasma- sheet  ions  of  ionospheric 
(rather  than  solar-wind)  origin  would  even  have  entered  the  plasma 
sheet  with  kinetic  energy  proportional  to  ionic  mass. 

It  is  not  really  accurate  to  say  that  the  observed  hardening  of 
the  proton  spectrum  with  decreasing  L  in  Figure  2  is  a  consequence 
of  Liouville1  s  theorem.  It  is  not  true,  for  example,  that  7(M.  J,L)  = 
7(M,  J.  L  ).  This  can  be  seen  from  (6).  On  the  other  hand,  Liouville  s 
theorem  is  not  violated  by  the  radial- diffusion  process.  The  distinc¬ 
tion  that  must  be  made  here  is  that  Liouville' s  theorem  recuires  the 


Figure  3.  Equatorial  pitch-angle  distributions  of  alpha  particles  (left 
panel)  and  protons  (right  panel)  in  the  same  four  energy  /nucleon  pass 
bands,  as  observed  on  spacecraft  OV1-19  (Blake  et  al.  ,  1973).  Sta¬ 
tistical  error  bars  (not  shown)  would  be  about  as  large  as  the  data- 
point  symbols  themselves. 


entity  that  diffuses  is  the  longitude- averaged  (i.  e.  ,  o-averaged  or 
drift-averaged)  phase- space  density,  whereas  the  longitude -dependent 
phase-space  density  remains  conserved.  The  dilemma  confronting 
the  space  researcher  is  that  f(M,  J,  I  .  ip)  is  difficult  to  observe  with 
high  resolution  because  particles  that  differ  even  slightly  in  M  or  J 
have  different  drift  periods  and  thus  become  separated  rather  quickly 
in  longitude  (<p).  The  primary  observable  manifestation  of  drift-phase 
organization  in  the  phase-space  density  is  the  drift-echo  phenomenon, 
and  even  this  is  discernible  for  only  a  few  drift  periods  following  an 
impulsive  change  in  the  magnetospheric  field  (e.  g.  ,  Lanzerotti  et  al.  , 
1971).  Nevertheless,  the  permanent  dispersal  of  particles  in  L  by 
such  a  Hamiltonian  process  definitely  implies  their  permanent  organiza¬ 
tion  in  ip  for  each  M  and  J  in  the  distribution,  and  such  drift-phase 
organization  really  is_  a  consequence  of  Liouville's  theorem  since  the 
drift  phase  is  canonically  conjugate  to  the  action  variable  (q/c)4>  for  a 
particle  of  charge  q  and  third  invariant  <b  (Roberts,  1967;  Schulz  and 
Lanzerotti,  1974). 


The  normalization  constants  7  X  10"^  day"*  and  1  X  10"*°  day"* 

that  appear  in  (3a)  and  (3b),  respectively,  are  those  obtained  by  Cro- 

ley  et  al.  (1976)  from  a  variational  analysis  of  inner-zone  proton  data. 

It  follows  that  dJJ  5  if  M/Z  2  120  MeV/gauss  (i.e.,  if  E  2 

3 

40  Z/L  MeV)  for  ions  that  mirror  on  the  equator  (i.e.,  for  y  =  1). 
For  y  «  1  one  finds  2  only  if  E  2  180  Z/L^  MeV.  These 

results  imply,  for  example,  that  ring-current  ions  (for  which  E  ~  10- 
100  keV)  are  transported  mainly  by  electrostatic  impulses  (unsteady 
magnetospheric  convection).  Among  ions  for  which  E  2  1.2  Z/L  MeV 


( e ) 


scales  as 


[so  that  ~  ^0  i-n  (3b)],  the  diffusion  coefficient  D  T • 

2  2  2 

(Z/A)“  at  fixed  K“  (=  J"  /Sm^M)  and  fixed  M/A  (first  invariant  per 
nucleon).  This  is  the  scaling  to  which  Cornwall  (1972)  appealed  in 
his  argument  for  preferential  access  of  H+  (Z/A  =  1)  over  He  (Z/A 
=  1/2)  and  He  (Z/A  =  1/4)  to  the  inner  part  of  the  magnetosphere, 
e.  g.  ,  to  L  <  3.  Conversely,  among  ions  for  which  E  <  30Z/L  keV 
[so  that  (f2^r)“  S  1/40  in  (3b)]  the  diffusion  coefficient  becomes  inde- 
pendent  of  M,  J,  Z,  and  A.  In  this  case  one  obtains  D^^/  a  8  X 
10  ^(t/20  min)^L^  day  Ring-current  ions  occupy  the  middle  range 

(0.  1  £  ^3t  -  10)  of  the  parameter  between  these  two  extremes 

and  thus  defy  efforts  to  simplify  the  expression  for  Perhaps 

this  is  not  an  accident,  inasmuch  as  the  decay  time  r  of  an  impulse 
in  the  convection  electric  field  should  naturally  reflect  the  time  scale 
of  the  underlying  particle  motion  whereby  the  convection  electric  field 
might  be  discharged. 


Although  analyses  based  on  radial  diffusion  alone  are  quite  in¬ 
structive,  additional  transport  processes  must  be  invoked  if  one  is  to 
account  for  further  details  of  the  prevailing  ionic  phase-space  distri¬ 
butions.  For  example,  one  may  ask  where  (in  L)  the  radiation  intensity 
peaks  at  a  given  energy  E.  Differentiation  of  (6)  with  respect  to  L 

j*; 

with  (9a)  inserted  for  7(M,  J,  L  )  yields  (97/3L)j.  =  0  under  the  con¬ 
dition 


7(L:/L)7 
1  -  (I^/L)7 


E  1  d  In  En 

f  +  1  -  -  - - 

Eq  J  d  In  L 


(10) 


I 


If  the  boundary  spectrum  is  assumed  to  be  Maxwellian  (so  that  (  = 
then  for  particles  mirroring  at  the  equator  one  obtains 


1 


for 


as  the  location  of  peak  radiation  intensity,  since  Eq  =  (L""/L)'Eq' 
y  =  1  in  (5)  and  (9b).  The  factor  [1  -  (L^/L)^]^^  varies  only 
from  0.99  to  1.00  as  L  varies  from  about  1.6  to  infinity,  and  so  the 
expectation  based  on  (6)  is  that  the  location  in  L  of  maximum  radia¬ 
tion  intensity  should  vary  as  E  *  .  Observations  compiled  by 

Fritz  and  Spjeidvik  (1981)  show  clearly  that  the  actual  variation  closely 
-3/16 

resembles  E  ,  as  had  been  predicted  by  Tverskoy  (1965),  and  not 

E  *  ^  for  1.5  £  L  £  4.  5;  moreover,  the  peak  alpha-particle  inten¬ 
sity  occurs  about  1.3  times  as  far  from  the  geocenter  as  the  peak 
intensity  for  protons  of  the  same  energy,  and  this  too  had  been  pre¬ 
dicted  by  Tverskoy  (1965).  The  essential  process  considered  by  Tver¬ 
skoy  (1965)  in  addition  to  radial  diffusion  was  Coulomb  drag,  which 
formally  causes  a  particle's  M  to  decrease  while  K  (  =  J  /8mpM) 
remains  constant,  i.  e.  ,  causes  a  particle's  energy  to  decrease  while 
the  equatorial  pitch  angle  remains  constant.  In  order  to  take  advan¬ 
tage  of  the  fact  the  K  is  conserved  by  both  processes  (radial  diffusion 
and  Coulomb  drag)  it  is  usual  to  transform  from  the  canonical  variables 
(M,  J,  9)  to  the  "new”  variables  (M,  K,  L).  The  Jacobian  of  this  trans¬ 
formation  has  an  absolute  value  given  by 

|d(M,  J.«)/d(M.K,L)|  =  (8m0M)1/2(2iT(i/L2a)  (12) 

and  so  the  appropriate  Fokker  -  Planck  equation  is  of  the  form 


Here  the  transport  coefficient  (dM/dt)  is  given  by 

(dM/dt)y  =  M‘1/2(4ffZ2q^ /me)(L9a9y6m0/2|j.3)1/2 

X<£N.Z.[y2  -  1  -  Y2ln(4iameM/m0y2L3a3I.)] 

+  Ne[Y2  -  1  -  Y2ln(\Dmev/1i)]  >  (14) 

where  Y^  =  1  +  (2|j.M/mQC2y2L3a3)  =  [1  -  (v/c)2]’*,  and  where  the 
subscript  v  denotes  "frictional  drag"  as  distinguished  from  other  pro¬ 
cesses.  Other  symbols  appearing  in  (14)  include  the  charge  qg  and 

rest  mass  m  of  an  electron,  the  densities  N  of  free  electrons  and 
e  e 

Nj  of  neutral  atoms  or  molecules  of  atomic  number  Z.,  the  mean 

energy  I.  required  for  atomic  excitation  or  ionization,  the  Debye 

length  X.£j,  and  the  symbol  "h  (Planck's  constant  divided  by  2ir).  Thus, 

for  energetic  ions  of  the  same  M/A  and  K  ,  the  transport  coefficient 

-1/22 

(dM/dt)v  scales  as  A  Z  .  This  means  that  the  ratio  M/A  decreases 
at  a  rate  proportional  to  A  3^2Z2  for  an  ion  of  charge  state  Z  and 
atomic-mass  number  A,  i.  e. ,  as  64:32:8:1  for  H+:  He++:  He+:  0+. 

The  scaling  of  (dM/dt)y  with  L  naturally  depends  on  the  models 

adopted  for  atmospheric  and  plasmaspheric  density  distributions.  Tver- 

skoy  (1965)  described  Coulomb  drag  in  terms  of  a  "lifetime”  propor- 
3/2  -2  -1/2  - 1 

tional  to  E  Z  A  throughout  the  radiation  belts,  and 

thereby  replaced  the  basic  transport  equation  with  an  equation  that 
would  yield  some  simple  results  in  closed  form.  The  model  atmosphere 
and  plasmasphere  described  by  Farley  and  Walt  (1971)  yielded  a  "life- 
time"  (  -  d  In  E/dt)  roughly  proportional  to  L  *  for  protons  of 
specified  M  and  vanishing  K  in  the  interval  1.  6  S  L  S  2.  5  (and  pre¬ 
sumably  for  higher  L  values  within  the  plasmaspher e ).  This  would 


correspond  to  a  transport  coefficient  (dM/dt)  ,  proportional  to  L“ 
over  a  region  in  which  hydrogen  atoms  and  free  electrons  are  the  dom¬ 
inant  agents  of  Coulomb  drag.  The  ion  energy  lost  to  neutral  atoms 
and  molecules  in  (14)  results  in  elevation  of  bound  electrons  to  higher 
energy  levels  (excitation)  or  to  the  spectral  continuum  (ionization). 

The  energy  lost  to  free  electrons  (via  the  term  containing  and  X^) 
results  in  the  Cherenkov  emission  of  Langmuir  waves.  The  logarithms 
in  (14)  are  typically  quite  large  (  a  10)  for  radiation-belt  ions. 

The  angle  brackets  in  (14)  denote  the  guiding -center  average 
over  the  bounce  and  drift  motion.  Evaluation  of  the  bounce  average 
typically  entails  a  numerical  computation,  but  a  few  analytical  esti¬ 
mates  are  available  in  closed  form.  Cornwall  (1975)  suggests  that 
the  bounce -aver aged  atomic-hydrogen  density  scales  as  y  ^  times  the 

equatorial  density.  Smith  et  al.  (1976)  suggest  that  the  bounce- 

3  5 

averaged  density  of  hydrogen  atoms  scales  as  sec  ’  X.  ,  where  X-m 
is  the  mirror  latitude.  Schulz  (1977)  has  described  a  method  of  cal¬ 
culating,  for  any  positive  integer  n,  the  bounce -aver aged  value  of 
(B/Bq)11  along  a  dipolar  field  line  on  which  Bq  a  )  is  the  equa¬ 

torial  value  of  the  magnetic -field  intensity  B.  For  example,  the  case 
n  =  1  yields  the  result 

<(B/Bq)>  s  (  1.  9191  y"5  74  -  1 .  1786  y" 1  ) 

4  (  1.  3802  -  0.6397y3/4)  (15) 

which  might  be  relevant  to  a  model  in  which  is  proportional  to  B 
along  a  plasmaspheric  field  line.  For  magnetic  drift  shells  that  in¬ 
tersect  the  model  plasmapause,  account  must  be  taken  of  the  relative 
amounts  of  time  spent  inside  and  outside  the  plasmasphere  in  evalu- 


ating  the  azimuthal  (drift)  averages  of  and  Ngln( X^m^v /■&)  for  use 
in  (14).  The  dependence  of  X^  on  is  not  a  severe  complication 
here,  since  the  argument  of  the  logarithm  is  typically  very  large,  i.  e 
since  the  dependence  of  the  logarithm  on  is  typically  quite  weak. 
Thus,  it  should  be  sufficient  to  define  Xq  as  the  value  of  X^  at  B/B 
=  <(B/Bq)>  and  to  approximate  the  bounce  average  of  Nfiln( X^n^v/li) 
by  the  quantity  <Ne>ln(  X^m^v/li). 


0 
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LOSS  PROCESSES 


If  Coulomb  drag  is  considered  a  transport  process,  as  seems 
appropriate,  then  the  main  loss  mechanism  for  ring -current  and 
raaia  "ion-belt  ions  appears  to  be  charge  exchange.  In  its  simplest 
manifestation  charge  exchange  causes  a  singly  charged  ion  to  become 
neutral  and  thus  escape  magnetic  confinement  by  the  earth's  field. 

This  process  is  described  formally  by  adding  a  term  -  7/r^  to  the 
right-hand  side  of  (13).  The  charge-exchange  lifetime  is  given  by 
=  <N^><rv,  where  is  the  density  of  atomic  hydrogen  in  the 
exosphere,  v  is  the  speed  of  the  ion  [  as  in  (14)  above],  and  or  is  the 
ion's  cross  section  for  charge  exchange  with  atomic  hydrogen.  The 
structure  of  (13)  makes  it  desirable  to  specify  as  a  function  of  K, 
L,  and  M.  Results  for  H  and  He  ions  at  K  =  0  in  a  model  atmos- 
phere  of  atomic  hydrogen  are  shown  in  Figure  4  (Cornwall,  1971). 
Lifetimes  for  H+  ions  are  thus  typically  shorter  than  for  He+  ions  of 
the  same  K  and  M/A.  Both  and  o-v  vary  with  L  so  as  to  pro¬ 

duce  the  results  shown  in  Figure  4;  the  ion  energy  (on  which  o-v  de- 

_3 

pends  quite  strongly)  varies  as  L  for  K  =  0.  The  products  erv  are 
shown  as  a  function  of  ion  energy  for  various  ions  in  Figure  5  (Tins¬ 
ley,  1976;  Cornwall  and  Schulz,  1979).  The  cross  sections  <r  corre¬ 
spond  to  neutralization  for  H+,  C+,  N+,  and  0+.  Those  for  helium 
ions  describe  this  as  well  as  other  possibilities,  viz.,  He+T  — ►  He 
and  He'  — ►  He^  (dashed  curves)  and  He+  — ►  He  +  (shown  by  the  J-J-- — 
curve  at  £  2  150  keV  in  Figure  5). 

The  possibility  of  two  or  more  charge  states  for  an  ionic  species 


Figure  4.  Charge-exchange  lifetimes  against  neutralization  in  model 
atomic -hydro gen  atmosphere  for  equatorially  mirroring  (K  =  0)  protons 
(H+.  left  panel)  and  helium  ions  (He+,  right  panel)  at  selected  values 
of  M/A,  measured  in.  MeV/gauss -nucleon  (Cornwall,  1971;  Schulz  and 
Lanzerotti,  1974). 


ENERGY,  keV 


Figure  5,  Products  of  ion  velocity  and  cross  section  for  charge  ex¬ 
change  in  atmosphere  of  atomic  hydrogen  (Cornwall  and  Schulz,  1979), 
baaed  on  compilations  summarized  by  Cornwall  (1971)  and  Tinsley 
(1976).  Cross  sections  correspond  to  neutralization  for  H+,  C+,  N+ 
(solid  curves),  0+  (dotted  curve),  and  He+  (dashed  curve);  to  reduction 
of  charge  state  (He++  to  He+)  for  He++  (dashed  curve);  and  to  increase 
of  charge  state  (He+  to  He++)  for  He+  (curve  of  "plus"  signs  at  lower 


requires  that  coupled  transport  equations  be  introduced.  For  helium 
ions  the  relevant  phase-space  densities  are  7^  and  7,,  the  phase-space 
densities  of  He+  and  He++,  respectively.  The  coupled  transport  equa¬ 
tions  are  of  the  form 

Oij/at)  +  CjTj  =  DjTj  -  (^/tjq)  -  (7j  /  t12)  +  ( 72  /t2  j)  (16a) 

and 

072/3t)  +  C272  =  D272  -  (?2/t21)  +  (7j/t12)  (16b) 

where  C  and  D  are  the  Coulomb-drag  and  radial-diffusion  operators 
that  appear  in  (13),  while  r.^  denotes  the  lifetime  against  transition 
from  charge  state  i  to  charge  state  j.  Thus,  the  lifetime  in  (16a) 

is  identical  with  the  that  appears  in  the  helium  (He+)  panel  of  Fig¬ 

ure  4.  Spjeldvik  and  Fritz  (1978)  have  generalized  (16)  for  application 
to  the  charge  states  of  oxygen,  which  are  quite  numerous  but  present 
no  difficulty  in  principle. 

Strictly  speaking,  the  source  terms  72/t21  in  (16a)  and  7j/tj2 
in  (16b)  should  be  replaced  by  nonlocal  source  terms  in  which  T.(L') 
is  averaged  (with  appropriate  weighting)  over  a  narrow  interval  j  L' 

-  L  |  <  AL  about  the  L  value  of  interest.  Such  a  refinement  of  (16) 
would  acknowledge  the  fact  that  a  change  in  charge  state  inherently 
displaces  the  guiding  center  of  a  geomagnetically  trapped  ion  by  in¬ 
stantaneously  altering  the  gyro-radius.  For  ions  of  initial  charge  state 
Z  and  initial  gyro- radius  p  at  the  equator,  the  maximum  change  in  L 
value  from  transition  to  charge  state  Z±1  would  be  of  order  AL  = 
(p/a)(Z±l)”*,  where  &  is  the  radius  of  the  earth.  Such  considerations 
would  normally  be  unimportant,  since  the  adiabatic  description  of 


charged-particle  morion  (upon  which  ring-current  and  radiation-belt 
theory  is  founded)  requires  p/La  «  y  £  1  for  its  validity.  However, 
charge  exchange  from  state  i  to  state  j  does  not  remove  an  ion  from 
the  radiation  belt  unless  j  =  0,  nor  does  charge  exchange  entail  sig¬ 
nificant  pitch-angle  diffusion  of  the  trapped  ion.  Thus,  the  same  ion 
can  experience  charge  exchange  between  the  same  pair  of  states  many 
times  during  its  residence  in  the  magnetosphere,  and  the  effect  of  in¬ 
troducing  nonlocal  source  terms  to  account  for  this  in  the  equation  for 
37./8t  would  simulate  augmentation  of  the  radial -diffusion  coefficient 
by  an  amount  AD./'P  ~  <  ( p/a)^( Z  ±  1  )~^(r. .  +  t. .)**>.  Such  an  augmen- 

1}  J1 

tation  might  (or  might  not)  be  quantitatively  significant. 

Transition  from  charge  state  Z  =  1  to  charge  state  Z  -  1  (=0) 
displaces  the  guiding  center  to  infinity  and  thus  violates  the  adiabatic 
description  of  charged-particle  motion.  Such  an  event  (i.  e. ,  neutrali¬ 
zation)  removes  a  fast  ion  from  the  ring  current  but  creates  a  fast 
neutral  whose  direction  of  motion  depends  on  a  random  variable,  viz.  , 
the  ion's  gyrophase  at  the  instant  of  neutralization.  It  takes  time  (~ 

20  sec)  for  the  fast  neutral  to  escape  from  the  magnetosphere,  and 
during  this  time  the  fast  neutral  can  be  "stripped"  of  one  or  more  e- 
lectrons  so  as  to  re-appear  as  a  born-again  ion  elsewhere  in  the  mag¬ 
netosphere.  Such  reconversion  of  a  fast  neutral  into  an  ion  is  unlikely 
to  occur  unless  the  trajectory  of  the  neutral  carries  it  through  the 
dense  atmosphere  immediately  surrounding  the  earth.  However,  Mi- 
zera  and  Blake  (1973)  have  indeed  reported  evidence  of  this  process 
in  the  form  of  a  partial  belt  of  born-again  ions  of  ring-current  energy 
at  L  <  1.  1  near  the  magnetic  equator.  Observation  of  the  effect  is 


facilitated  by  the  eccentricity  (~  440  km)  of  the  drift  shells  relative 
to  the  earth,  since  ions  reborn  at  a  100-km  altitude  just  west  of  the 
South  Atlantic  anomaly  would  gradient-drift  to  an  altitude  >  900  km 
in  about  half  a  drift  period. 

Despite  the  above-described  implications  for  ion  transport,  charge 
exchange  is  primarily  a  loss  mechanism  for  ring -current  ions.  It  fol¬ 
lows  from  Figure  5  that  H+  is  lost  more  rapidly  than  0+  or  He+  at 
ion  energies  <  50  keV.  Lyons  and  Evans  (1976)  have  emphasized  the 
importance  of  charge  exchange  in  this  context  for  the  ion- composition 
of  the  recovery -phase  ring  current.  Their  conclusion  was  that  the 
ring-current  ions  observed  on  the  satellite  (also  known  as  Explorer 
45)  decayed  much  too  slowly  in  intensity  (and  developed  pitch-angle 
anisotropy  much  too  slowly)  to  have  been  H+  ions.  This  was  a  dis- 
turbing  suggestion  at  the  time,  since  the  predominance  of  H  among 
the  ions  of  the  ring  current  had  by  then  become  an  article  of  faith. 

The  argument  of  Lyons  and  Evans  (1976)  was  not  air-tight,  since  it 
ignored  (a)  the  likely  possibility  that  radial  diffusion  would  continue 
during  recovery  phase  and  (b)  the  adiabatic  energization  that  ring- 
current  particles  would  experience  through  the  azimuthal  electric  field 
induced  by  decay  of  the  ring  current  itself.  However,  the  work  of 
Lyons  and  Evans  (1976)  alerted  space  researchers  to  the  likelihood 
(now  generally  accepted  as  fact)  that  ions  heavier  than  H+  are  major 
constituents  of  the  ring  current  at  least  part  of  the  time.  The  abun¬ 
dance  of  0+  in  auroral  ion  beams  (see  above  under  "Access")  makes 
this  conclusion  seem  less  surprising  in  retrospect. 

The  creation  of  pitch- angle  anisotropy  in  ring-current  ion  aistri- 


butions  is  an  expected  consequence  of  charge  exchange,  since  the 
bounce-averaged  density  of  atomic  hydrogen  in  the  atmosphere 

should  be  an  increasing  function  of  mirror  latitude  at  each  L 

value.  It  was  noted  above  in  connection  with  Coulomb  drag  that  sev¬ 
eral  functional  forms  have  been  proposed  for  scaling  with  equa¬ 

torial  pitch  angle  (o^  m  sin  *y)  or  mirror  latitude,  e.g.,  y  *  by 

Cornwall  (1975)  and  sec^'^X  bv  Smith  et  al.  (1976).  The  similaritv 

m  - - •— 

of  these  two  scalings  is  revealed  by  the  relationship 

y  =  (1  +  3  sin^X.)**^  cos^X  (17) 

that  holds  between  y  and  X^  in  a  dipolar  field.  If  processes  other 
than  charge  exchange  are  neglected  for  simplicity,  then  an  initially 
isotropic  phase-space  density  7(E,L;0)  evolves  into  the  distribution 

T(E,  L,  y;t)  =  exp(  -  <NH>«rvt)  I(E,  L;0)  (18) 

which  is  anisotropic  for  times  t  >  0  by  virtue  of  the  dependence  of 

on  y.  Cornwall  (1977)  realized  that  such  anisotropy  among  ring- 

current  ions  could  lead  to  instability  in  electromagnetic  ion-cyclotron 

wave  modes.  He  was  able  to  calculate  the  growth  rate  in  closed  form 

as  a  function  of  wave  frequency  (w/2tt)  and  elapsed  time  (t)  under  the 

simplifying  (but  admittedly  unrealistic)  approximation  that  scales 

-2  - 1  35 

as  y  rather  than  as  y  or  as  sec  '  X  .  His  results  are  neverthe- 
1  m 

less  revealing  and  qualitatively  consistent  with  expectation.  The  growth 
rate  is  initially  negative  at  all  frequencies  w/2t r  >  0  because  the  pitch- 
angle  distribution  is  isotropic,  as  can  be  seen  from  (18)  for  t  =  0. 
Thereafter  the  distribution  becomes  anisotropic,  and  so  the  growth  rate 


becomes  positive  for  a  band  of  frequencies  ui /2~  between  2ero  and 
some  fraction  of  the  ion  (in  this  case,  proton)  gyrofr equency  £2/ 2”. 
Experience  shows  that  the  fraction  should  be  approximately  A  /(A  -  1), 
where  A  is  the  anisotropy  of  the  distribution.  The  fraction  is  exactly 
A /(A  +  1),  for  example,  if  the  pitch-angle  distribution  has  the  form 
sin  ffg  or  corresponds  to  a  bi- Maxwellian  distribution  for  which  the 
temperature  ratio  is  given  by  T^/T1(  =  A  +  1.  Of  greater  interest 
than  the  local  growth  rate  I*  is  the  path-integr ated  gain 

G  =  exp  (2r/vg)  ds  (19) 

where  v^  is  the  group  velocity  and  ds  is  the  element  of  arc  length 
along  the  ray  path  (usually  taken  to  be  a  field  line  in  order  to  simplify 
the  task).  It  is  common  in  magnetospheric  physics  to  assume  that  G 
a  exp(4rQLa/Vg),  where  rQ  is  the  equatorial  growth  rate,  in  order  to 
simplify  the  task  further.  However,  Eiemohn  (1967)  has  demonstrated 
interesting  effects  that  are  overlooked  if  the  path  integral  is  not  evalu¬ 
ated  more  carefully.  In  any  event,  the  gain  G  has  a  maximum  with 
respect  to  w/2ir,  and  the  maximum  gain  thus  determined  is  not  mono¬ 
tonic  as  a  function  of  t.  Indeed,  the  value  of  Gm  (t)  is  larger  than 
unity  for  0  <  t  <  oo  but  equal  to  unity  for  t  =  0  (when  the  anisotropy  A 
vanishes)  and  for  t  *  oo  (when  I  itself  varnishes  by  virtue  of  charge  ex- 

change).  Thus,  if  G  (t)  increases  monotonically  for  t  <  t"  and  de- 

max 

creases  monotonically  for  t  >  t*,  the  remaining  question  is  whether 

G  (t*)  is  sufficiently  large  for  instability  to  occur  at  all.  The  usual 
max 

criterion  for  instability  is  G_  (t  )  >  400,  i.  e.  ,  a  single-pass  gain  of 

m  aX 

at  least  20  in  wave  intensity  (or  13  dB  in  decibel  notation)  from  one 
foot  of  the  field  line  to  the  other,  since  it  is  presumed  that  only  about 


3 tj  of  the  wave  intensity  that  reaches  either  foot  of  the  field  line  will 
be  reflected  back  into  the  magnetosphere.  The  integral  in  (19)  refers 
to  propagation  once  in  each  direction  along  the  field  line  and  so  must 
compensate  for  two  such  5%  reflections.  If  the  instability  occurs,  i.  e. 

3*e 

if  G _ (t  )  i  400,  then  the  resulting  ion- cvclotron  waves  will  enhance 

the  loss  rate  of  ring-current  ions  by  causing  pitch-angle  diffusion  into 
the  loss  cone.  Such  diffusion  will  naturally  reduce  the  anisotropy  of 
the  pitch- angle  distribution  and  thus  limit  the  instability,  but  it  will 
also  enhance  the  effectiveness  of  charge  exchange  as  a  loss  mechanism 
for  ring-current  ions  by  increasing  the  mirror  latitude  (and  hence  the 
magnitude  of  <N^>)  for  the  average  ion  trapped  in  the  flux  tube. 

It  follows  from  the  above  considerations  that  the  probability  of 
having  pitch- angle  diffusion  among  ring -current  ions  increases  with  the 
strength  of  the  magnetic  storm,  i.  e.  ,  with  the  intensity  of  the  ring 
current  itself.  Recovery  from  a  weak  storm  might  be  achieved  through 
charge  exchange  alone,  whereas  recovery  from  a  major  storm  might 
involve  both  charge  exchange  and  pitch-angle  diffusion  (each  serving  to 
enhance  the  effectiveness  of  the  other).  The  foregoing  scenario  is  in¬ 
complete  in  at  least  three  respects.  First,  the  initial  pitch- angle  dis¬ 
tribution  is  not  likely  to  be  isotropic,  since  the  radial  transport  that 
populates  the  ring  current  from  the  plasma  sheet  creates  pitch- angle 
anisotropy  even  when  charge  exchange  and  Coulomb  drag  are  neglected. 
This  was  illustrated  above  (but  for  higher -energy  ions)  in  Figure  2  and 
in  equations  (5)-(9).  Thus,  wave  growth  occurs  even  at  t  =  0,  although 

G _  (0)  is  likely  to  be  smaller  than  G  (t  )  for  some  t  >  0. 

max  max 

Secondly,  the  ring  current  is  likely  to  be  populated  by  two  or 


more  ionic  species  simultaneously.  Perhaps  K'  is  the  major  ion  ini¬ 
tially.  i.  e.  .  during  the  main  phase  of  a  magnetic  storm,  with  He”, 

He  ,  and  O  present  as  (relatively)  minor  constituents.  In  this  case 
the  electromagnetic  mode  structure  is  interrupted  by  stop-bands  above 
the  gyr ©frequencies  of  the  heavier  ions,  and  the  hydrogen-cyclotron  in¬ 
stability  is  partially  suppressed  in  favor  of  enhanced  excitation  at  fre¬ 
quencies  below  the  gyrofrequencies  of  the  helium  and  oxygen  ions  (e.  g. , 
Cornwall  and  Schulz,  1971;  Chiu  et  al.  ,  1980).  Ions  having  gyrofre- 
quency  £1/ 2tt  can  experience  pitch-angle  diffusion  through  a  Doppler- 
shifted  cyclotron  resonance  with  a  frequency  w/2ir  represented  in  the 
wave  spectrum  if  the  ion  energy  is  sufficient  to  satisfy  the  resonance 
condition 

kv  =  (u  -  fi)sec  a  >  0  (20) 

where  2ir/k  is  the  wavelength  that  corresponds  to  the  wave  frequency  w/2n 
Resonance  for  Q  >  w  requires  cos  a  <0,  i.  e.  ,  particle  and  wave  trav¬ 
eling  in  opposite  directions  along  the  field  line.  Resonance  for  D  <  w 
requires  particle  and  wave  to  travel  in  the  same  direction  along  the 
field  line  so  that  cos  a  >  0.  As  H+  ceases  to  be  the  major  ion,  the 
stop  bands  above  the  helium  gyrofrequencies  grow  larger  and  suppress 
the  hydrogen- cyclotron  instability  altogether  (e.  g.  ,  Chiu  et  al.  ,  1 9 . ) . 

J.+  4-  + 

but  instability  can  still  occur  in  the  He'  ,  He  ,  and  O  passbands. 

Finally,  the  foregoing  scenario  is  incomplete  in  that  electrostatic 
waves  have  been  neglected.  The  usual  rule-of-thumb  in  magnetospheric 
physics  is  that  electromagnetic  waves  predominate  inside  the  plasma- 
sphere  and  electrostatic  waves  outside.  Electromagnetic  ion- cyclotron 
waves  are  easier  to  analyze  because  the  cold-plasma  dispersion  rela¬ 
tion  can  be  perturbed  in  order  to  extract  hot-plasma  effects  such  as 


wave  growth.  This  does  not  usually  work  for  electrostatic  ion- 
cvclotron  waves,  since  in  this  case  the  cold-plasma  dispersion  relation 
does  not  provide  a  good  approximation  of  the  phase  velocity.  The 
usual  spectrum  of  elect? ostatic  ion-cyclotron  turbulence  is  characterized 
by  excitation  bands  located  about  midway  between  harmonics  of  the  ion 
gyrofreauency,  and  the  usual  source  of  free  energy  is  anisotropy  (e.  g.  , 


T^  >  T|()  in  the  pitch-angle  distribution.  Interesting  complications 
should  ensue  for  the  ring-current  plasma  outside  the  plasmasphere, 
since  (for  example)  the  j  harmonic  of  the  H  gyrofrequency  coincides 


PERSPECTIVE 


Implementation  of  the  foregoing  principles  may  impress  the  reader 
as  an  onerous  task.  It  is  natural  to  ask  whether  a  simpler  theory  of 
magnetosp'neric  ion  composition  is  not  already  available.  In  a  certain 
sense  the  theory  of  Tverskov  (1965)  fills  this  request.  Tverskoy  (1965) 
argued  that  the  maximum  in  radiation  intensity  at  any  energy  E  should 
be  located  at  the  L  value  for  which  the  "time  scales"  associated  with 
radial  diffusion  and  Coulomb  drag  are  equal.  This  would  occur  where 
M  *(dM/dtty  =  -  L  2D^^  in  the  language  of  the  present  work.  Tver¬ 
skoy  (1965  )  neglected  the  Coulomb  drag  associated  with  neutral  atoms 
and  molecules  in  (14),  as  well  as  the  radial  diffusion  associated  with 
electrc  itatic  impulses  (unsteady  convection)  in  (3).  On  this  basis  one 
would  expect  to  find 

(4trZ2q^/me)(m0/2E3)1/2<Ne>ln(kDmev/1i) 

«  8.  1  X  10  "  14  L8  [Q(y)/180  D(y)]2  sec"1  (21) 

at  the  L  value  of  maximum  radiation  intensity  for  specified  E,  Z,  A, 
and  y,  since  the  argument  of  the  logarithm  in  (21)  is  very  large  (> 

1 0 "^ )  and  the  value  of  y  is  of  order  unity  in  most  situations  of  interest. 
Solution  of  (21)  for  L  yields 

Lmax  *  3.00<(Ne/103  cm'3)>1/S  [(l/16)ln(XDmev/-{j)]1/8 

X  [180D(y)/Q(y)]1  74  Z1 /4  A1 /16  (E/1  MeV)'3716  (22) 

as  the  location  of  peak  radiation  intensity.  The  normalization  of  N  , 

_  3  -3  . 

E,  and  ln(  X^m^v/li )  by  typical  values  (10  cm  ,  1  MeV,  and  It,  re- 
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spectively)  of  these  quantities  makes  the  above  expression  for  L 

easier  to  evaluate  by  inspection.  It  was  on  the  basis  of  a  similar 
expression  that  Tverskoy  (1965)  anticipated  the  E"  scaling  law  for 

Lffiax  and  the  appearance  of  the  alpha-particle  (Z  =  2,  A  =  4)  maximum 

3  /8 

at  1.3  (»  2  )  times  as  high  an  L  value  as  the  intensity  maximum 

for  protons  of  the  same  energy. 

The  scaling  laws  deduced  by  Tverskoy  (1965)  are  amply  confirmed 
by  the  observations  that  Fritz  and  Spjeldvik  (1981)  have  recently  com¬ 
piled.  Even  the  leading  factor  3.00  in  (22),  which  is  not  arbitrary 
but  derived  from  (21),  seems  appropriate  in  that  (for  y  =  1)  the  1-MeV 
proton  maximum  occurs  at  L  -  3  and  the  1-MeV  alpha-particle  maxi¬ 
mum  at  L  4  (actually  at  L  =  2.76  and  L  =  3.59,  respectively,  upon 
more  careful  inspection).  Such  good  agreement  seems  surprising  in 
view  of  the  fact  that  (21)  and  (22)  have  no  fundamental  basis  but  only 
heuristic  justification.  On  closer  examination,  however,  the  agreement 
is  less  impressive  since  (for  example)  the  E  scaling  law  for 

follows  from  (21)  only  if  is  independent  of  L.  The  more  likely 

alternative  scaling  in  which  <Ng>  is  proportional  to  L  would  have 
yielded  at  Z  *  ^  A*  E~  *  in  substantially  poorer  agreement 

with  the  data  compiled  by  Fritz  and  Spjeldvik  (1981).  Furthermore, 
the  "alpha-particle"  observations  ordinarily  *  **fer  to  a  mixture  of  Z  =  1 
and  Z  s  2  charge  states  rather  than  to  Z  =  2  alone,  since  the  charge- 
exchange  rates  that  prevail  in  the  earth's  magnetosphere  make  this 
mixture  unavoidable  (Cornwall,  1972).  It  is  not  the  purpose  of  this 
discussion  to  quibble  over  details  of  significant  work  that  is  nearly 
twenty  years  old,  but  rather  to  illustrate  how  much  more  sophisticated 
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the  treatment  of  radiation-belt  dynamics  has  become  over  the  intervening 
period.  Just  as  the  heuristic  estimates  by  Tverskoy  (1965)  satisfied  the 
needs  of  an  earlier  decade  of  space  research,  so  the  work  of  Spjeldvik 
and  Fritz  (1978)  exemplifies  the  more  fundamental  and  quantitative  ap¬ 
proach  that  is  demanded  today. 

The  present  work  is  not  intended  as  a  compendium  of  ''who  has 
done  what"  in  the  field  of  magnetospheric  ion  composition,  but  rather 
as  an  introduction  to  the  principles  upon  which  the  modern  dynamical 
theory  of  geomagnetic  ally  trapped  ion  populations  is  based.  Review 
articles  with  extensive  bibliographies  on  the  subject  have  been  prepared 
by  Cornwall  and  Schulz  (1979),  Spjeldvik  (1979),  and  Fritz  and  Spjeld¬ 
vik  (1981),  among  others.  The  present  assignment  seemed  to  require 
a  different  approach,  and  the  foregoing  pages  represent  the  outcome 


of  this  endeavor. 
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